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O ■ Abstract 
(N 

^ ■ The amplitude for the anomalous process 77r — >■ tttt is evaluated to two loops in the 



ix. 



chiral expansion by means of a dispersive method. The two new coupling constants 
that enter at this order are estimated via sum rules derived from a non-perturbative 
chiral approach. With these coupling constants fixed, the numerical results are given 
and compared with the available experimental information. 
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1 Introduction 



At low energies, the strong interaction between pions can be described by the 
^ ■ effective field theory of QCD called chiral perturbation theory (ChPT) [1-4]. 

. In this effective field theory the interaction between pions is analyzed in terms 

of a systematic expansion in powers of the low external momenta and the 
small pion mass. This chiral expansion is equivalent to an expansion in the 
number of loops and works very well for processes involving pions [5]. 

In the normal intrinsic parity sector, the chiral expansion starts at O(p^) and 
is obtained from the leading order chiral Lagrangian. In this sector the one- 
loop or O(p^) corrections were extensively treated in the works by Gasser and 
Leutwyler [2,3], and many processes have now been calculated to this order. 
The further extension to two loops of O(p^) has also been undertaken with 
several two-fiavour calculations already published in the normal sector [6-10]. 

For the abnormal intrinsic parity sector, the chiral expansion starts at 0(p'^) 
and is obtained from the Wess-Zumino-Witten effective action [11]. In this 
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sector the one-loop corrections of 0{p^) have also been analyzed [12], and 
some anomalous processes have been calculated to this order [13]. However, 
contrary to the case in the normal sector, no two-loop calculations of 0{p^) 
have been published for processes in the anomalous sector. 

In this paper the anomalous process 771 tttt is calculated to two loops using a 
dispersive method. This process is important for the theory of chiral anomalies 
and has previously been calculated to both leading [14] and one-loop [15] order 
in the chiral expansion. However, both of these results are somewhat below 
the present experimental data [16]. This has motivated new and more precise 
experiments, which will be performed at different facilities such as CERN [17], 
FNAL [18], and CEBAF [19]. 

It is therefore important to evaluate the two-loop corrections and compare the 
result with the experimental information. This is the purpose of the present 
paper, which is organized as follows. In Sect. 2, the notation and kinematics of 
the process — * tttt are given, whereas previous ChPT results are reviewed 
in Sect. 3. The general dispersive formula is derived in Sect. 4, together with 
the explicit calculation of the two-loop corrections. In Sect. 5, this two-loop 
result is used in a non-perturbative chiral approach in order to estimate the 
two new coupling constants from sum rules. The numerical results are given 
in Sect. 6, including the comparison with experiments, and Sect. 7 consists of 
a short conclusion. 



2 Kinematics 



The amplitude for the anomalous process 

7(g)7r"(pi) 7r"(p2)7r°(po) 



is given in terms of the scalar function ^3^(5, t, u) as 



(2) 



where s = {pi + q^, t = {pi — ^2)^, and u — (pi — PoY are the Mandelstam 
variables. In the following it is assumed that all particles are real. In this case 
s + t + u — 3M^, so that in the center of mass system one has 
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3M^-s+{s-M^)^l- AMl/s cos 6 
?>Ml-s-{s- M^)^l- AM^/s COS0 



(3) 
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with the center of mass scattering angle. Because of isospin symmetry, the 
other 77r — > tttt reactions are given by the same scalar fmiction ^3,^(5, t, -u), 
which is fully symmetrical in its arguments. The total cross section is obtained 
from the expression 

^^'^ = 1024^ ,1/2 (^-^^') I desiT.H\F,^{s,t,u)\\ (4) 



and the partial wave expansion of the scalar function F^t^{s, t, u) is of the form 

F3,(s,i,ii)= ^/Ks)P/(cose). (5) 

odd I 

The projection onto the lowest P partial wave is given by 

h{s) = ^J dQsm'eFs^{s,t,u), (6) 

and the higher partial waves can be projected out in similar ways. This partial 
wave expansion contains only odd I since the pions in the final state have the 
total isospin / = 1. In the unitarity relation, the complete set of intermediate 
states are the same as for tttt scattering. Within the elastic approximation 
containing only the two-pion intermediate state, the unitarity relation is 

lmMs)^a{s)fns)tns) (7) 



with a{s) = y {s — 4:M^)/s the phase-space factor and t} the nn partial wave 
with isospin I — 1 and angular momentum I. This unitarity relation implies 
that the phase of /; will coincide with the tttt phase shift Sf in accordance 
with Watson's final-state theorem [20] . Inelasticities due to other intermediate 
states such as the four-pion or the KK states remain very small below 1 GeV 
and are completely negligible at low energies. 



3 Chiral expansion 

Anomalous processes start at O(p^) in the chiral expansion and are obtained 
from the Wess-Zumino-Witten effective action at tree level. With only the 
electromagnetic interaction as external fields, these anomalous processes con- 
tain exclusively an odd number of pseudoscalars, i.e., they have abnormal or 
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odd intrinsic parity. For the process 777 — > vrvr, the leading order chiral result 
for the scalar function F2T^{s,t,u) is given by [11,14] 

F3.(5,t,«) = Fi? = ^|^, (8) 



where A^c = 3 is the number of colors and F^^ — 92.4 MeV is the pion decay 
constant. Another related anomalous process is the decay tt^ —>■ 77, where the 
Ft^jj coupling obtained at leading order is of the form [11,21] 

ir(0) ^ ^ (Q) 



These two results are exact in the soft pion and soft photon limit where they 
are related to each other by the low-energy theorem = F^^^/eF^. For the 
decay 7r° — > 77, the prediction from Eq. (9) is in excellent agreement with the 
experimental information, thus providing important evidence for the three- 
color nature of the strong interaction. However, for the process 77r tttt, 
the corresponding prediction Fa^r = 9.7 GeV~^ is somewhat lower than the 
experimental measurement Fa^ = 12.9 ± 0.9 ± 0.5 GeV~^ [16], which could 
even be said to favor the value Nc — 4^. 

It is therefore of importance to calculate the corrections to the soft pion and 
soft photon limit. Within ChPT these corrections can indeed be calculated in a 
systematic manner. The first correction is the one-loop contribution of O(p^), 
which has previously been calculated for both processes. In the case of the 
decay 7r° 77, these one-loop corrections turn out to be very small [22], i.e., 
they do not spoil the excellent agreement with the experimental information. 
On the other hand, for the process 77r — > tttt, the one-loop contributions are 
larger and they will therefore be of importance when comparing with the 
experimental measurement. The expression for the scalar function F^T^[s,t,u) 
to this order in the chiral expansion is [15] 



(0) 



647r2 

—-Cl{,,){s + t + u) + C]! 



(10) 



where the term Cf contains the contributions from the loops and is given by 



C7 ^ 



967r2F2 



Ml 

-{s + t + u)\og^ + -{s + t + u) 
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with F{x) expressed in terms of the standard one-loop function J{x) as 



F{x) = 4:7r^{x - 4:Ml)J{x) - \x. 



(12) 



The expression (10) converges to the chiral anomaly in the soft pion 
and soft photon limit and is given in terms of the renormalized low-energy 
constant from the 0{p^) anomalous chiral Lagrangian [12]. This low-energy 
constant depends on the renormalization scale ^ and is needed to absorb 
divergences in the one-loop calculation. In principle, C2 should be determined 
phenomenologically, preferably from other observablcs, but at present this 
appears to be rather out of reach. Therefore, this low-energy constant has been 
estimated using the assumption of vector resonance saturation, which is known 
to work well for the 0{p^) non-anomalous low-energy constants [23]. Assuming 
that the same is the case for the 0{p^) anomalous low-energy constants, one 
has [15] 



with the renormalization scale typical chosen at the resonance scale ji — Mp — 
770 MeV. Having fixed the value of this low-energy constant, it is possible 
to obtain the prediction for the one-loop expression (10). This improves the 
agreement with the experimental measurement [16]. However, the prediction 
is still on the lower side of the data [15,17,18]. 

Therefore, it is important with both theoretical and experimental improve- 
ments. Indeed, significant improvements on the experimental side are expected 
from several new experiments [17-19]. A similar improvement on the theore- 
tical side would involve the calculation of the two-loop corrections. This could 
be done by a full field theory calculation to two loops. However, these two-loop 
corrections can also be calculated by a dispersive method, which will be the 
subject of the next section. 



4 Dispersive representation 

J,..! Derivation of the dispersive formula 

The anomalous process 77? nn is in many ways similar to the tt^tt^ 7r°7r° 
scattering process. They are both described in terms of a single scalar function, 
which is fully symmetrical in the s, t, and u variables. The tttt scattering 
process can be described by Roy equations [24] derived from the fundamental 
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1287r2M2 



5 



principles of analyticity, crossing, and unitarity. When these Roy equations are 
combined with the chiral expansion, the general structure of the irn scattering 
amplitude can be obtained from a dispersive representation to two loops in 
the chiral expansion [8,25]. 

The same is also the case for the process 771 tttt. In order to show this, 
one starts with a fixed-t dispersion relation for the scalar function ^3,^(5, t, u). 
Fixed-t dispersion relations with two subtractions are known to exist for tttt 
scattering due to the Proissart bound. Assuming that the same is also the case 
for the process 77r — > tttt, one has the fixed-t dispersion relation 

F3,(., t) = C(t) + i / i^'M^fA + ^) , (14) 

TT J S'^ \S' — S S' — U I 

4M2 ^ ^ 



where u = 3M^ — s — t. The s ^ u symmetry of ^3,^(5, t, u) implies that there 
is no subtraction constant linear in s. Therefore, the only subtraction constant 
is C(t), where the t dependence can be obtained from the s ^ t symmetry 
-^37r(0,t) = F3^(t, 0). This implies that C{t) can be written as 



00 

C(i)=C(0) + - / ds' 



,ImF3,(s',0) f ^ ml-ty 



s'2 \s'-t s'-^Ml + t 



00 



1 7 ^^ MF^^js'.t) mi-ty .... 

^ I 7^ s'-3M2 + r ^^^^ 

AMI 

which can be inserted in Eq. (14). By construction, the resulting dispersion 
relation exhibits s <-> symmetry for fixed whereas s t symmetry for 
fixed u is not manifest. In order to impose this latter symmetry one can expand 
the absorptive part of F3^(s',t) in partial waves using Eq. (5). Writing this 
expansion in the form 

ImF3,(s', t) = Im/i(s') + M(s', t), (16) 



the higher partial waves with I > 3 are contained in the function Im$(s',i), 
which is given by 

Im$(s',t) = J2^mfi{s')P;{cose). (17) 

l>3 

With this decomposition of the partial waves, it is possible to rewrite Eq. (14) 
with C{t) given by Eq. (15) as 

F37r(s, t, U) = F3^{S, t, U) + $37r('S, t, m), (18) 
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where ^3^(5,^,14) and ^3T^{s,t,u) are given by 



F,As, t, u) = m + - / ds'^ 7^ + 7^ + 7^ ' 

4M2 ^ ^ 

4M2 ^ / 

This is the Roy equation for the anomalous process 77r — > tttt. One observes 
that -F37r(s, t, m) is now fully symmetrical in s, t, and u, whereas this symmetry 
is not manifest in $37r(s, t, u). However, at low energies, the absorptive part of 
the higher partial waves with Z > 3 is negligible. This implies that in practice 
^z-nis-it^u) can be treated in a simple way as a small, real correction at low 
energies. This fact makes the corresponding Roy equations for tttt scattering 
very useful. 

In ChPT the absorptive part of the higher partial waves is indeed suppressed. 
Within this methodology, these higher partial waves only start at 0{p^) in 
the chiral expansion. Since the corresponding tttt partial waves start at O(p^), 
perturbative unitarity implies that the absorptive part of the higher partial 
waves with / > 3 is of 0(p^°) or higher in the chiral expansion. Thus, up to 
this order, the term $37r(s,t, m) in Eq. (18) vanishes, and the full amplitude 
is given entirely by the function F^t^{s, t, u). In order for the full amplitude to 
formally satisfy the chiral anomaly in the soft pion and soft photon limit, one 
may write the subtraction constant as C(0) = -F3°^[l + C{s + t + u)], where 
s + 1 + u = 3M^. Thus, assuming for the moment that two subtractions are 
indeed sufficient, the full amplitude to 0{p^) in the chiral expansion may be 
written as 



(0) 



+ 



TT J 

4M2 



l + C{s + t + u) 



ds' 



u 



-,12 



S — S 



s' — t s' — U 



(20) 



In this dispersive representation, the absorptive part of the lowest partial wave 
can be determined from unitarity. In the general unitarity relation, the 2n-pion 

invariant phase space is of 0(p^""^), the amplitude for multi-pion scattering 
is dominantly of O(p^), and the amplitude for multi-pion photo-production is 
at least of 0{p*). Consequently, intermediate states containing more than two 
pions are suppressed at least up to 0{p^^) in the chiral expansion. Therefore, 



7 



within SU(2) ChPT, the absorptive part of the lowest partial wave is given by 
elastic unitarity to 0{p^) as 



where /{" is the P partial wave of 0(p^"+^) and ti the corresponding tttt 
partial wave of 0(p^"+^) [2]. In SU(3) ChPT, one must also include the in- 
elasticity from the KK intermediate state. However, since this inelasticity is 
completely negligible at low energies, only SU(2) ChPT will be considered in 
the following. 

4-2 The amplitude to two loops 

From the dispersive representation (20), one obtains straightforwardly the one- 
loop ChPT formula (10) by setting Im/i = lmf[^\ In this case the subtraction 
constant C can be expressed in terms of the low-energy constant C2 and chiral 
logarithms. However, in order to calculate the amplitude to two loops, it is 
necessary to use one subtraction more in the dispersive representation. This 
is due to the fact that the absorptive part of fi^^ behaves as modulo log 
factors. Thus, the amplitude to two loops can be obtained from the following 
dispersion relation: 



where C and D are the two subtraction constants and Im/i is given by 



Im/i = Im/i^^ -|- Im/i^-' with Im/}^-' and Im/}^-' determined from Eq. (21). 



This dispersion relation can be evaluated with the same methodology as has 
been applied in the calculation of the tttt scattering amplitude [8] and the pion 
form factors [6] to two loops. With the use of this methodology, the result can 
be written as 



lm/r(«) = 0, 

Imfi'\s)^a{s)fi'\s)tl^'\s), 




(21) 



F3^(s, t, u) = F^l> [l + C{s + t + u) + D{s^ + t^ + u^) 
1 7 , Mf lis') ( 




(22) 



F3^{s,t,u) 



F^'J [l + Cis + t + u) + D{s^ + e + u^) 

+t/^(s) + U^{t) + [/^(«) + [/^(s) + ^7^(t) + U^{u)\ , (23) 



where the U'^ term can be expressed in the compact analytic form: 
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(24) 



The explicit expressions for the functions Ki{x) and K/i{x) are given in Ref. [8], 
where these functions were introduced in the evaluation of the tttt scattering 
amplitude to two loops. They are analytic functions with cuts starting at the 
TTTT threshold and can be expressed in terms of the standard one-loop function 
J{x). As for the other term U^, this part does not have a compact analytic 
representation in terms of elementary functions. However, it can be obtained 
numerically from the dispersive representation: 



U^{x) 



X 



IT 



ds' 



AM2 



s'^{s' — x) 



(25) 



In this dispersive formula t^i^^ is the lowest order tttt P partial wave and 
is given by 

^^^'^ = hl ^96^^ {^^^ + + ^t^^^^ + ^^""A ^^^^ 



with t and u determined from Eq. (3) and F{x) given by Eq. (12). The size of 
the corrections from are very small at low energies and they can therefore 
in practice be neglected at these energies. 

In the dispersive representation, Eq. (23), the full amplitude to two loops 
is determined up to the subtraction constants C and D. These subtraction 
constants may be parameterized in terms of the coupling constants ci, C2, and 
o?2 as 



C- 



167r2F2 



IX Ml _' 
(^-1 - e) + ^^^^ 
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D- 



1 



1 



+ 



60M2 167r2F2 



(27) 



where the couphng constant ci can be expressed in terms of the 0{p^) anoma- 
lous low-energy constant C2 and chiral logarithms: 



(28) 



The couphng constants C2 and ^2 enter at two-loop order in the chiral ex- 
pansion and contain contributions both from two-loop diagrams and from the 
unknown renormalized low-energy constants, which parameterize the O(p^) 
anomalous chiral Lagrangian. 



5 Non-perturbative chiral approach 



5. 1 The lowest partial wave 



Prom the two- loop expression for the amplitude, Eq. (23), one may project 
out the lowest P partial wave with the use of Eq. (6). This gives the expansion 

h{s) = f?\s) + f['\s) + J?\s), (29) 



which will satisfy the pcrturbativc unitarity relations (21). These relations 
work very well at low energies, whereas the deviation from exact unitarity, Eq. 
(7), becomes more pronounced as the energy is increased. At these energies, 
still higher order unitarity corrections will start to be of importance in the 
chiral expansion. This is particularly the case in the p(770) resonance region, 
where unitarity corrections are essential. 

In the presence of the p(770) resonance, unitarity will therefore be of the 
utmost importance. There are different ways to combine exact unitarity and 
the chiral expansion in order to try to account for this resonance. One such 
method, which has been successfully applied to many different processes, is 
the so-called non-perturbative inverse amplitude method (1AM) [26-29]. Since 
this is a rather general method, it can also be straightforwardly applied to the 
present case. The starting point for the JAM is to write down a dispersion 
relation for the inverse of the partial wave /i. In this dispersion relation, exact 
unitarity and the chiral expansion are used to compute the important right 
cut, whereas the left cut and the subtraction constants are approximated by 
ChPT. If one writes down a similar dispersion relation for the chiral expansion 
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using perturbative unitarity on the right cut, it is possible to express the result 
of the lAM in a simple way in terms of the chiral partial waves. With the two- 
loop ChPT expansion (29), the result of the lAM can be written as 



Ms) 




(30) 



This is formally equivalent to the [0,2] Pade approximant applied to ChPT 
and will therefore coincide with the chiral expansion up to two loops. However, 
since exact unitarity was used in the derivation of the lAM, it is expected that 
this result will improve ChPT at higher energies. The 1AM applied to two- 
loop ChPT has been extensively discussed in Ref. [28], where the detailed 
derivation of the general result can also be found. 

The result of the lAM, Eq. (30), depends on the pion mass and pion decay 
constant, which will be set equal to Mj^ = 139.6 MeV and = 92.4 McV, 
respectively. Furthermore, the 1AM is also given in terms of the low-energy 
constants I2 — h and I4. These low-energy constants appear in the ChPT tttt P 
partial wave to one loop order and can been determined phenomenologically 
from other sources. However, since the 1AM contains higher order unitarity 
corrections, the phenomenological values of these low-energy constants in the 
1AM do not necessarily coincide precisely with the values obtained in ChPT. 
For the combination I2 — li, this has been determined by applying the lAM 
to one-loop ChPT in the case of tttt scattering with the result ^2 — = 5.8 
[28]. The other low-energy constant I4 has been determined in ChPT from the 
pion scalar form factor with the central value I4 = 4.4 [7]. Since the 1AM (30) 
does not depend much on the precise value of this low-energy constant, the 
same value can also be apphed in the lAM. Finally, this result also depends 
on the one- loop couphng constant Ci. This coupling constant is related to the 
low-energy constant C2, which has been estimated from the assumption of 
vector resonance saturation, Eq. (13) [15]. With the value of chosen at the 
renormalization scale jj, = Mp, one finds Ci = 1.71. Therefore, in the lAM (30), 
the values 



will be used throughout. Having fixed these constants, it is possible to deter- 
mine the remaining two-loop coupling constants C2 and ^2 in the 1AM from 
the mass and width of the p(770) resonance. The mass of this resonance can 
be defined as the energy where the phase passed 90° and the width can be 
determined from the slope of the phase at the resonance. With Mp — 770.0 
MeV and = 150.7 MeV [30], this gives the values 



l2-h^ 5.8 



h = 4.4 



ci = 1.71 



(31) 



C2 = -27 



d2 = 3.477, 



(32) 
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Fig. 1. The phase of the lowest partial wave /i obtained from the lAM, Eq. (30), 
compared to the experimental vrvr phase shift Si from Ref. [31] (circles), Ref. [32] 
(squares), and Ref. [33] (diamonds). 




Vs (GeV) 

Fig. 2. The lowest partial wave obtained from the lAM, Eq. (30). 

where the lAM depends rather strongly on d2 and to a lesser extent on C2- 
With these values the lAM gives the phase shown in Fig. 1. It is observed that 
the result agrees very well with the experimental phase shifts all the way up to 
1 GeV. Therefore, the lAM satisfies unitarity at least up to this energy. In Fig. 
2, the normalized absolute square of the JAM partial wave is shown. In this 
case there is no experimental data to be compared with. However, assuming 
that the full amplitude is completely dominated by the lowest partial wave in 
the resonance region, the cross section (4) can be expressed in terms of 
This cross section is usually parameterized in the resonance region in terms 
of the p — > 7r7 partial width F^^j by using the Breit-Wigner form 
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Equating the two expressions for the cross section at the resonance energy, 
the value F^r^ = 96 KeV is obtained from the lAM. Experimentally, the value 
of has been extracted from the production of p via the Primakoff effect by 
incident pion on nuclear targets. The present experimental data are, however, 
not very consistent with each other [30]. The latest experiment gave the value 
= 81 ± 4 ± 4 KeV [34], whereas two earlier experiments gave somewhat 
lower values for F^r^ [35]. Prom the process e'^e~ — > 777, this partial width has 
also been obtained giving the value F^,, = 121 ± 31 KeV [36]. Therefore, at 
present it can only be concluded that the value for F^^ obtained from the JAM 
is not in conflict with the experimental situation, but new experiments in the 
resonance region are definitely needed. 



5.2 The full amplitude 

If the absorptive part of the higher partial waves with Z > 3 is neglected, it 
was shown in Sect. 4.1 that the full amphtude could be constructed from the 
absorptive part of the lowest partial wave. Since the lAM gives a good de- 
scription of this partial wave up to a least 1 GeV, one may use the absorptive 
part of this partial wave in the Roy equation for 77r tttt, Eq. (20). With the 
lowest partial wave approximated by the lAM, the Roy equation may be writ- 
ten with only one subtraction, which is fully determined by the chiral anomaly 
result. Therefore, the non-perturbative chiral result for the full amplitude can 
be written as 

F3,(.. «) = Fi? + - J + ' + , (34) 

n J s' \s' — s s' — t s' — uj 

4M2 

where Im/i is given by the absorptive part of the lAM result (30). The full 
amplitude is by construction fully symmetrical in s, t, and u. However, if the 
lowest partial wave /i is projected out from Eq. (34), the result does not agree 
exactly with the lAM result (30). This is due to the fact that the left cut is 
now fully determined from crossing symmetry instead of being approximated 
by the chiral expansion. In the elastic region below 1 GeV, the difference is, 
however, negligible. From the full amplitude, Eq. (34), it is also found that the 
higher partial waves are very small below 1 GeV. Therefore, the cross section 
can indeed be expressed in terms of to a very good approximation. For a 
somewhat different evaluation of the Roy equation for 777 —> tttt, see Ref. [37]. 
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As the lAM was derived using elastic unitarity, it is expected that this result 
is only applicable in the clastic region. The high energy part of the dispersion 
relation (34) is therefore not expected to be very well approximated by the 
lAM. Still, at low and moderate energies the dispersion relation should be 
almost completely saturated by the p(770) contribution. Therefore, at these 
energies the high energy part of the dispersion relation should not be very 
important. Furthermore, at these energies it is also expected that the contri- 
bution from the higher partial waves to the Roy equation should be negligible. 
Therefore, the non-perturbative chiral result (34) should give a rather accu- 
rate description of the full amplitude at low and moderate energies. Prom this 
result, it is possible to determine the subtraction constants C and D in the 
two-loop ChPT result (23) via the sum rules 

1 1 7 ^. wi(.o ^_ 1 1 7 ^ /^h{s') .35. 

37r 4jv^2 Stt 4jv^2 

Evaluating these sum rules, one obtains the values C = 0.93 GeV~^ and 
i) = 2.12 GeV~^, respectively. These sum rules are indeed almost completely 
saturated by the dispersion relation up to 1 GeV, as the higher energy part 
only contributes by approximately 2% for C and approximately 1% for D. 
With the one-loop coupling constant Ci determined from the assumption of 
vector resonance saturation, Eq. (31), this gives the following values for the 
two-loop coupling constants 

C2 = -20 , J2 = 2.7. (36) 

These values are slightly different from the values determined directly from 
the lAM (32). As it has already been stated, this is indeed to be expected 
since the lAM contains higher order unitarity effects, which will effect the 
determination of these coupling constants. 

Another way to combine unitarity and the chiral expansion has been proposed 
by Holstein [38]. In this model, vector meson dominance is combined with 
one-loop ChPT in order to account for rescattering effects. Projecting out the 
lowest partial wave, it is found that unitarity is approximately satisfied and 
that the p — > 7r7 partial width is given by P^-y = 84 KeV. Within this model, 
one obtains the values 

C2 = -19 , (i2 = 2.9, (37) 

which agree very well with the determination given in (36). However, from the 
model proposed by Holstein, all higher partial waves have the same phase as 
the lowest partial wave, which is in contradiction with the experimental facts. 
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In order to remedy this, the lowest partial wave obtained in this model may 
be inserted in the dispersion relation (34). The resulting amplitude still agrees 
well with the experimental information, but now implies that all the higher 
partial waves are real. This gives the values C2 = —16 and d2 = 2.8 from the 
sum rules (35), which is also fully consistent with the values given in (36). 
Therefore, it seems that the solution of the sum rules is rather robust against 
small variations of the input. 

In the numerical results for two- loop ChPT (23), the values of C2 and J2 
determined in (36) will be used together with the value of Ci given in (31). For 
the low-energy constants li and I2, these have been determined in ChPT. With 
the use of the two-loop result for the tttt scattering amplitude and fitting the 
D partial wave scattering lengths, the value h — h = 6.0 has been obtained [9]. 
Similar values have also recently been obtained from a two-loop calculation 
of the Ki4 form factors [39] . These values are rather consistent with the value 
given in (31) but are somewhat lower than the value h — h = 7.8 obtained 
from a dispersive improved one-loop calculation of the K14 form factors [40]. 
However, the numerical result only depends very slightly on the exact value 
of I2 — h- Therefore, the values given in Eq. (31) for both I2 — h and U will 
also be used for ChPT. 



6 Numerical results 

6.1 Amplitude and cross section 

Fig. 3 shows the absolute square of the amplitude normalized to the chiral 
anomaly. In the sub-threshold region, the two-loop result is rather close to 
the one-loop result, whereas the corrections start to be of importance slightly 
above threshold. As for the one-loop result, this gives an almost constant cor- 
rection to the leading order chiral anomaly result. This is due to the fact that 
the only variation comes from the s, t, and u dependency of the function F 
in the one-loop expression (10). The contribution from the subtraction con- 
stants is also shown in the figure, where this contribution is obtained from 
the expression F3^{s,t,u) = Fj°^[l + C{s + t + u) + D{s^ + + m^)]. It is 
observed that this part gives the main contribution to the two-loop result, 
which is due to the presence of the p(770) resonance. However, the unitarity 
corrections coming from the U(x) terms are also of some importance in the 
numerical result. In fact, these unitarity corrections are essential in the non- 
perturbative chiral result, Eq. (34), which is also shown in the figure. The 
deviation between this result and the two-loop chiral result gives an estimate 
of the range of validity of the truncated chiral expansion. From this deviation, 
it is observed that still higher order chiral corrections should begin to be of 
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2 4 6 8 10 12 14 16 

Fig. 3. The amplitude \F^t^/F^\'^ evaluated at t = u as a function of s/M'^. The 
solid line is the two-loop result, Eq. (23), the dashed line the one-loop result, Eq. 
(10), and the dotted line the leading order chiral anomaly result, Eq. (8). Finally, the 
long-dashed line is the non-perturbative chiral result, Eq. (34), and the dash-dotted 
line is the contribution from the subtraction constants. 

some importance somewhat around s — lOM^, which is also observed in other 
two-loop calculations. 

In Fig. 4, the total cross section is shown as a function of s/M^, where the total 
cross section is given by Eq. (4). Since the one-loop expression for F3^(s, t,^) 
is more or less constant in the whole phase space, the one-loop result for 
the cross section gives an almost constant correction to the leading order 
result of approximately 20 %. Close to threshold, the additional two-loop 
correction is rather small compared to the one-loop correction, whereas this is 
not the case for higher energies. The additional two-loop correction amounts 
to approximately 75 % of the one-loop correction at s = lOM.^, and even more 
at higher energies. This indicates that even higher order terms in the chiral 
expansion should begin to be of importance at this energy. This is also observed 
in the figure, where the non-perturbative chiral result begins to deviate from 
the two-loop result around this energy. 

6.2 Comparison with experiments 

The process 777 — > tttt has been investigated at low energies by the Serpukhov 
experiment [16]. This experiment used the Primakoff reaction of pion pair 
production by pions in the nuclear Coulomb field 

TT- + (Z, A) ^ TT- + 7r° + (Z, A) . (38) 
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Fig. 4. The total cross section as a function of s/M^. The solid line is the two- loop 
result, the dashed line the one- loop result, the dotted line the leading order result 
and the long-dashed line the non-perturbativc chiral result. 

The cross section for this process is related to the 77r — > tttt cross section 
through the equivalent photon method 



da 



dsdtdq^ 



TT 



1 da. 



77r— >7r7r 



s-Ml dt 



(39) 



where 



dt 



\F3^{s,t,u) 
5127r 



(40) 



and 



2E 



(41) 



and E is the energy of the incident pion beam. Neglecting the dependency 
in F3T,{s,t,u), i.e., setting ft! so that s + t + u — 3M^, the total cross 
section is given by 



a — 



ds 



:s-4M2)3/2 



10247r2 J sV2 
dO sire' 9\F^T,{s,t,u 



I 



In 



2 

?max 



X 



(42) 
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Table 1 

The total cross section for the Primakoff reaction (38) divided by Z^. The theo- 
retical predictions are obtained from Eq. (42) with F3^(s,t, u) given by the leading 
order chiral anomaly result {0{p^)), the one-loop result (0(p^)), the two-loop result 
{0{p^)), and the non-perturbative chiral result (NPCR). The experimental result is 
the average value obtained in Ref. [16], where the first (second) error is statistical 
(systematic). 



0(p4) 


0(p6) 


0(/) 


NPCR Experiment 


a/Z^ (nb) 0.92 


1.09 


1.18 


1.21 1.63 ± 0.23 ± 0.13 



Table 2 

The chiral anomaly extracted from the experimental data. The different va- 
lues correspond to setting F^.^(s,t,u) in Eq. (42) equal to the leading order chiral 
anomaly result (0(p*^)), the one-loop result {0{p^)), the two-loop result {0{p^)), 
and the non-perturbative chiral result (NPCR). The theoretical result is from Eq. 
(8) with = 92.4 ± 0.3 MeV [30]. 



0(p4) 


0(/) 


0(/) 


NPCR 


Theory 


^3^°^ (GeV-3) 12.9 ± 1.4 


11.9 ±1.3 


11.4 ± 1.3 


11.3 ± 1.3 


9.7 ±0.1 



The Serpukhov experiment was carried out with a 40 GeV pion beam and the 
maximum momentum transfer was r/f^.^^ = 2 x lO^'' GeV^ <^ M^. Thus, the 
dependency in ^3^(5, t, u) can indeed be neglected when one compares with 
this experiment. 



The kinematical region studied in the Serpukhov experiment was restricted 
to Sinax = lOM^. Three different targets (C, Al, and Fe) were used and the 
measured total cross sections were found to agree well with the theoretical Z'^ 
dependence. Averaging the values of a/Z'^ obtained from the three different 
targets, the experimental result given in Table 1 was obtained. This experi- 
mental result has to be compared with the theoretical predictions, which are 
also given in Table 1. It is observed that the two-loop contribution gives a 
sizable correction to the one-loop result, whereas the further correction from 
the non-perturbative chiral result is rather small. Nevertheless, all results are 
still somewhat below the central experimental value, which, however, also has 
rather large error bars. 

Prom the experimental result of cr/Z^, it is possible to extract the value of the 
chiral anomaly ^ if this is regarded as a free parameter. These extracted 
values are shown in Table 2, where the statistical and systematic errors have 
been added. The theoretical prediction obtained from Eq. (8) is also shown in 
this table. With ^37^(5,^,^) given by the leading order chiral anomaly result, 
one obtains a value of F^ which is 2.3(T too high compared with the theoretical 
prediction. This result for F^ is the value generally quoted from the analysis 
in Ref. [16]. However, with the one-loop expression for ^3^(5, t, m), one obtains 
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a value of F^J which is only 1.7a too high, and with the two-loop result the 
disagreement is only at the 1.3a level. 

Since the non-perturbative chiral result is close to the two-loop result, the 
uncertainty coming from yet higher orders in the chiral expansion should be 
rather small, at least in the kinematical region probed by the Serpukhov ex- 
periment. There is of course also some uncertainty due to the uncertainty in 
the determination of the couphng constants C2 and ^2. However, as aheady 
discussed, this uncertainty should also be rather small. Therefore, in view 
of the disagreement between the theoretical predictions and the Serpukhov 
experiment, new improved Primakoff experiments are definitely needed. 

The process 771 tttt has also been investigated at low energies at CERN 
[41]. This experiment used the reaction 7r~e — > 7r~7r°e with 300 GeV pions and 
measured the total cross section for this process. It was found that the cross 
section obtained agreed with the chiral anomaly prediction with three number 
of colors. However, since the error bars were rather large, it was not possible 
to observe any systematic deviations from the soft pion and soft photon limit 
in this experiment. 

Therefore, new precision experiments are necessary in order to investigate the 
process 77r — > tttt and thus the theory of chiral anomalies in greater detail. 
Indeed, such new experiments are under way at different facilities. In the 
COMPASS experiment at CERN [17] and in the SELEX experiment at FNAL 
(E781) [18], the Primakoff reaction (38) will be measured with 600 and 50-280 
GeV pion beams, respectively. In these two experiments, the expected number 
of near threshold two-pion events is several orders of magnitude higher than 
previously obtained. This will allow analysis of the data separately in different 
intervals of s with small statistical errors. 

The SELEX experiment will also measure the reaction n^e ^ 7r~7r°e in or- 
der to obtain independent information on the 77r — > tttt amplitude. For this 
reaction, the expected number of events is also significantly larger than pre- 
viously obtained, which would give excellent complementary information on 

Finally, at CEBAF the process jti — > Tin is investigated by measuring 7p — > 
7r"'"7r''n cross sections using tagged photons [19]. Since the resonance region 
will also be measured in this experiment, this could give new information on 
the p — > 7r7 partial width F^-j,. However, as the incident pion is virtual in the 
CEBAF experiment, this has to be taken into account when comparing with 
theory. 
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7 Conclusion 



The anomalous process 77r — > tttt plays an important role in the theory of 
chiral anomalies. At leading order in the chiral expansion, which corresponds 
to the soft pion and soft photon limit, the amphtude for this process is given 

in terms of the number of colors Nc of the strong interaction. Comparing the 
leading order result with the present experimental information [16], one finds 
that the value A^c = 4 is favored. 

In order to test this important conclusion more precisely, it is necessary with 
both experimental and theoretical improvements. Indeed, significant improve- 
ments on the experimental side are expected from several new experiments 
[17-19]. On the theoretical side, the one-loop correction to the leading order 
result has previously been calculated [15]. However, in view of the new pre- 
cision experiments, it is important also to calculate the additional two-loop 
correction to the amplitude. 

This has been the purpose of the present paper, where the amplitude for the 
anomalous process ■jtt — > tttt has been evaluated to two loops in the chiral 
expansion by means of a dispersive method. The two new coupling constants 
that enter at two-loop order were determined from sum rules with the use of 
a non-perturbative chiral approach. The uncertainty in the numerical results 
due to this determination was estimated to be rather small. Moreover, the still 
higher order terms in the chiral expansion were also estimated to be small at 
low energies. 

The two-loop result improves the agreement with the present experimental 
information [16] compared to both the leading order and the one-loop re- 
sults. However, the two-loop prediction is still significantly below the central 
experimental data for Nc = 3. This fact is not likely to be due to theoreti- 
cal uncertainties. Therefore, should the new experiments [17-19] confirm the 
present central experimental value with better precision, it would be a serious 
problem for QCD. 
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